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1 Introduction: 



The Whitham equation 

/+oo 
K(x - y)u y dy = (1) 
-00 

was introduced in [I] to model a wave equation containing both breaking and peaking. A good 
discussion of equations ([I]) can be found in [2j. Here we are interested in equation (P) with the 
kernel 

K(x-y) = -\x-y\ , 

which yields 

ut + uu x + d^u = , (2) 

or in a local form 

(u t + uu x ) x + u = . (3) 

Equation ([3]) also follows as a particular limit of the following generalized Korteweg-de Vries (KdV) 
equation 

(ut + uu x - (3u xxx ) x = 7U , (4) 

derived by Ostrovsky [3] to model small-amplitude long waves in a rotating fluid (71* is induced 
by the Coriolis force) of finite depth. For (3 = (no high-frequency dispersion) the equation 
([3]) is known under different names in the literature, such as the reduced Ostrovsky equation, the 
Ostrovsky-Hunter equation, the short-wave equation and the Vakhnenko equation. From now on we 
will call (|3|) the Ostrovsky- Vakhnenko (OV) equation. This OV equation describes the short-wave 
perturbation in a relaxing medium [4]. This equation has a purely dispersive term and although it 
has the same nonlinearity of the KdV equation the dispersive terms are different. In [4] and in a 
series of papers El [7] it was established its integrability by deriving explicit solutions. Also, in 
[8] the integrability via inverse scattering method was derived via a third-order eigenvalue problem 
obtained after a Backlund transformation. However in [9] this Lax pair was written in its original 
variables as a zero curvature condition 

Ai,i - A 0rE - [A , Ai] = , 
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Also, they obtained the following third order Lax pair for the OV equation 



+ 



l 



u xx + 



1 



9A ' 3/ ° ' 

tpt ~ 9Xip xx + u^ x - u x ^ = , 

whose compatibility condition yields the x-derivative of the OV equation ([3]). 



(5) 



It turns out that the OV equation can be obtained as a short-wave limit of another integrable 
equation, the Degasperis-Procesi (DP) equation [lOl HH E2] 



u t ~ u xx t + Auu x = 3u x u xx + uu x 



(6) 



The authors of |lH[T2j found a Lax pair, derived two infinite sequences of conserved charges for ([6]) 
and proposed a bi-Hamiltonian formulation. Hone and Wang [9] have shown and explored the fact 
that the OV equation can be obtained as a limit of the DP equation through the transformation 



T : < 
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in the short wave limit e — > 0, 



DP OV - e 2 (u t +4uu x 



Using this short wave limit Hone and Wang obtained the scalar linear problem © from the corre- 
sponding Lax pair of the DP equation. Also, from the bi-Hamiltonian study of the DP equation 
performed in [11] they proposed two Hamiltonian operators for the OV equation but not in the 
original variable u. 

In this paper our main interest is to investigate the Hamiltonian integrability for the OV equa- 
tion directly in its original variable and evolutionary nonlocal form ([2]). As far as we known this is 
the simplest equation involving nonlinearity an nonlocality and will provide a good "laboratory" 
for study of nonlocal equation. 



Results on the OV equation frequently rely on the well known reciprocal transformation to the 
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special case of Ito's equation which can be written in Hirota's bilinear form [6j. We point out 
that a transformation between the OV equation and the Bullough-Dodd-Tzitzeica equation is also 
known [13] but not explored in the literature. In this paper we will introduce and explore a new 
third transformation. In Section [2] we find a fifth order symmetry to the OV equation and we show 
in Section [3] that the evolution equation associated to this symmetry can be transformed to the 
Caudrey-Dodd-Gibbon-Sawada-Kotera equation (CDGSK) [HI [151 EE] via a chain of Miura-type 
transformations. In Section [2] we also find the first Hamiltonian structure for the OV hierarchy. 
The transformations introduced in Section [3] transform dependent as well as independent variables 
and we give in Section 0] a formula relating the Hamiltonian structures in this situation. This 
formula is used in Section [5] to obtain the second Hamiltonian structure of the OV equation via the 
well known CDGSK's Hamiltonian structure. Our results are purely algebraic and no analytical 
justifications are provided. Questions about wave breaking, blow-up rates, boundary conditions, 
well-posedness and so on can be found, for instance, in [T7] and references within. 



2 Symmetries and a first Hamiltonian structure 

First let us look for some symmetries of the OV equation (|2|). It is well known \18\ [19] that a 
symmetry a of the evolution equation ut = K(u) should satisfy 



where 



and 



[a,K](u) = a'(u)[K(u)} - K'(u)[a(u)] 



(7) 



F'{u)[X] = jF(u + eX) 



6=0 



(8) 



is the Frechet or directional derivative of F(u) in the direction of the vector field X. In the case 
that F(u, v) depends on two variables we write F u and F v for the Frechet derivative with respect 
to u and v, respectively. If H is a functional then 



5H 

H'(u)[X] = (gr a dH,X) = (—,X 



(9) 



where (•, •} denotes the usual scalar product between the dual spaces. In what follows, we indicate 
by "*" the transpose of a operator with respect to the duality. From ([7]) it is easy to check that 



(T = , (Ti = K(u) , a 2 = u x , 03 



(1 + 3u xx ) 



-2/3 



(10) 
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are symmetries for the autonomous OV equation ([2]). This means that the OV equation and the 
local nonlinear equation 



3„, / 1 



-2/3 



ut = -^D 3 x \ [Uxx + -j . (11) 

where the right-hand side is the generator 03 of the higher symmetries and where the coefficient 
—3/2 is taken for simplicity, are members of the same hierarchy of equations and share many 
properties such as integrability. As far as we know this Ostrovsky-Vakhnenko fifth-order equation 
(OVF) (jlip is a new integrable equation but we will show that it can be transformed to the CDGSK 
equation in the next Section. 

By construction the OV equation ([2]) and the OVF equation (jlip share the same conserved 
charges as well Hamiltonian structures. Let us find the first Hamiltonian structure for these equa- 
tions |2Uj . Introducing the Clebsch potential u = <f> x the equation ([2]) can be written as 

4>xt + (t>x4>xx + <i> = . 

This equation can be obtained from a variational principle, 5 J dtdx C, with the Lagrangian density 

C = \Mx + \<g - \<P 2 ■ (12) 

This is a first order Lagrangian density and we can use, for example, Dirac's theory of constraints 
to obtain the Hamiltonian and the Hamiltonian operator associated with (I12p . The Lagrangian is 
degenerate and the primary constraint is obtained to be 

$ = 7T - -<f> x , (13) 

where tt = dC/d<f>t is the canonical momentum. The total Hamiltonian can be written as 
H T = j dx (7T0i - C + A$) = J dx ~<j>l + ^4> 2 + A (tt - hj) x \ 



(14) 



where A is a Lagrange multiplier field. Using the canonical Poisson bracket relation 

{0(x),vr(y)} = 5(x-y) , (15) 

with all others vanishing, it follows that the requirement of the primary constraint to be stationary 
under time evolution, {<3?(x), H?} = 0, determines the Lagrange multiplier field A in (j!4j) and the 
system has no further constraints. Using the canonical Poisson bracket relations (|15|) . we can now 
calculate 



K(r. ,/) = {$(x), <%)} = l -D y 5{y - x) - ^DJ(x - y) . (1(5) 
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This shows that the constraint (|13p is second class and that the Dirac bracket between the basic 
variables has the form 

{<Kx),<l>(y)} D = {4>(x),<f>(y)} - [ dzdz'{4>{x),<S>{z)}J{z,z'){®{z'),<l>{y)} = J(x,y) , 



where J is the inverse of the Poisson bracket of the constraint (|16|) , 



dz K(x, z)J(z, y) = 8(x — y) 



This last relation determines D x J(x,y) = —5{x — y) or J(x,y) = T>5{x — y) where 



V = -D 



-i 



(17) 



We can now set the constraint (fT3|) strongly to zero in (fb^) to obtain 



Using (|17p and the transformation properties of Hamiltonian operators (see (|37|) ). we get 

V = D x (V) (D x )* = D x , 

and the OV equation ([2]) can be written in the Hamiltonian form as 

8Hx 



u t = V v 



Su 



Hi 



dx 



(18) 



It can be easily checked that H\ is conserved by both the OV and OVF equations for rapidly 
decreasing or periodic boundary conditions. Taking the second x derivative of the OV equation we 
obtain the trivial conserved charge 



dxu xx , (19) 



and from this equation times (u xx + 1/3) 2 / 3 we also get 

9 f ( n 1/3 
H 2 = -- / dx I u xx + - 

as conserved charge for both the OV and OVF equations. In this way the OVF equation has the 
following Hamiltonian representation 



u t = V ^-^- ' Vx = D x , H 2 = ~ f dx (u xx + i 



1/3 
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Multiplying the OV equation ([2]) by u and integrating we also get the following conserved charge 

/ dxu 2 . (20) 



3 Transformation to the CDGSK equation 



The link between the KdV and modified KdV equations through a Miura transformation is not an 
isolated result in the theory of integrable models. Very often a nonlinear equation is equivalent to a 
known and well studied equation via a chain of Miura-type transformations and now we will show 
that we can relate the OVF equation with the CDGSK equation using this procedure. We follow 
the methods used in [211 [221 1251124] . 

The transformation 

(x, t, u{x, t)) ^ (x, t, v{x, t)): v = {u xx + 1/3)- 1 / 3 (21) 
relates the OVF equation with 

vt = ~v*D 5 x v 2 . (22) 

The right-hand side of the transformation (|2ip is obtained from the separant of the OV equation. 
(The separant of an evolution equation u t = f(u,u x , . . . , Ur n \) is df /dur n \, where ur n ) is the highest 
derivative of u with respect to x.) Let us note that (|22p is one of the non-constant separant evolution 
Fujimito-Watanabe equations (see |21j and references within). The equation (j22h has a separant 
v 5 and by the Ibragimov substitution 

(y,t,w(y,t)) * (M,«GM)) = { ^J^J)^ ( 23 ) 

is related with the constant separant equation 

w t = w 5y - 5Wy 1 w 2y W4y + 5wy 2 wl y w 3y , (24) 

where Wk y = dyW, k = 2,3,4,5 (and similar notations for derivatives are used in what follows). 
Finally, the transformation 

(y, t, w(y, t)) ^ (y, t, z(x, t)) : z = -Wy l w 3y (25) 

relates with the CDGSK equation 

Zt = z 5y + 5zz 3y + 5z y z 2y + 5z 2 z y . (26) 
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In summary, we have schematically 



B(u,v) = B(w,v) = B(w,z) = 

OVF — > (22) < — (El — > CDGSK, 



where the implicit transformations are 

B(w, z){y) = (z + w^wsy) (y) , 
B(w, v)(y) = (v(w) - w y ) (y) , 

B(u, v)(x) =(v- (u xx + 1/3)- 1 / 3 ) (x) . (27) 

4 Hamiltonian structure behavior under change of variables 

Transformations between two evolution equations generate transformations between the corre- 
sponding structures such as recursion operators, Hamiltonian structures, conserved charges and 
so on [25\ 126] . Let be the change of variables among dependent and independent variables 

f y = P(x,u (n h , 

, / (28) 
\ v = Q{x,u^) , 

where represents all derivatives of u with respect to x of order at most n. We want to relate 
the Hamiltonian representations 

Ut=V (u)° V (u) Eu{hh (2Q) 

OU 

~ hB ~ 

v t = V^°— = V^E v (h), (30) 
ov 



where 



H[u] = J dxh(x,u {n) ) , (31) 
H[v]= [ dyh(yM n) ), (32) 



and E u (/i) is the Euler operator acting on the Hamiltonian density h. The transformation (|28p 
defines an implicit function B(u,v) = 0, and we have 

B u u t + B v v t = , 

or 

v t = -Tut , (33) 
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where 

T = B- 1 B U , (34) 

and B u and B v are the Frechet derivatives defined in ([8]) . Now we use the following result (see [18] , 
Exercise 5.49, pg. 386) for the relation between the action of the Euler operator under a change of 
variables 

E u (h) = OE v (h) , (35) 

where 

0(R) = Q* U (D X P ■ R) - P:(D X Q ■ R) . (36) 

From (1331). (1291) and (1351) 

vt = -TV^OE v (h) 
and comparing with (|30p we finally get 

jy(v) = _ TV (u) _ ^ 

When the independent variable is not transformed B u = —Q u , B v = 1, T = —Q u and O = Q^. 
Also, when we calculate how a recursion operator R = T>2T>^ 1 transforms under f)28|) the operator O 
drops out. These are the usual results commonly found in the literature [251 [26]. The generalization 
of (|37p for a number of dependent and independent variables greater than one is straightforward. 

5 Second Hamiltonian structure 

We use the results of the last Section to obtain the second Hamiltonian structure of the OV equation 
([2]) from the known Hamiltonian structure of the CDGSK equation (j26H given by [16] 



D« = Dl + 2(zD y + D y z) , #« = J dy Qz 3 - . (38) 

From ([27]) and {55]) we have 

2% = Au; t , A = -w^ 1 ^ + Wy 2 w 3y D y = - V - l D x v 3 Dl , (39) 

v t = Bw t , B = D y — Wy 1 w 2y = v 2 D x v~ l , (40) 
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1 / 1\~ 4/3 1 

v t = Cut, C = --(u xx + -\ D 2 x = --v 4 D 2 x . (41) 

From (|39p ^ (|4ip we already see that is very convenient to use the variable v(x, t) because this 
factorizes all the operators involved, via the relations D y = vD x , w y = v, w 2y = vv x , w% y = v(vv x ) x , 
z = —(vv x )x and u xx = v~ s — 1/3; note also that 

D 2 y + z = D x v*D x v- 1 , 

D\ + AzD y + 2z y = v^D^D^D^- 2 . (42) 



Now using (|35p and f)36|) we have 
6HM T <yfl"W 



A—j^— , A = A* = D^Wy 1 - D yWy 2 w 3y = v D 2 x v 3 D > x v- 2 , (43) 

SEW - n 

— = B— — , B = -D yWy - w 2y = -D x v z , (44) 

5h1> <?i, C-C—U(^ + \)^—\l&. (45) 



5u Sv ' 3 x V 3 y 3 

Therefore, the Hamiltonian operators transform as 

£>H = A~ X V^A' X , pW = BV^B , 1?W = CT X V { ^C- 1 , (46) 

and as result of (02J), (f39]) - (|4T|) and (f4~3"j) - (}4"5j) we obtain the operator of order minus five 

= gz)- 2 ^- 2 ^^- 1 ^ 3 ^- 1 ^^- 2 ^- 2 , (47) 



where now t> is not a dependent variable but just a placeholder for the expression in (|2ip . i.e., 
t)E( Ura + l/3)- 1 /3, 

Now, let us find the Hamiltonian H^ u ' corresponding to H^ z ' in (|38p . Its variational derivative 

is 

5H^ 1 2 

and from (1431) and (1441) we can transform it into 



5HW = i ,(5irW 

— — = B A— — = 

ov oz 

= v 3 v 6x + 9v 2 v x v 5x + I7v 2 v 2x v ix + \Avv 2 x v ix + -^v 2 v\ x + 

20 5 

+34w a .« 2a; U3 a ; - 2ti 3 u 3x + — - 10vlv% x + hv~ v x v 2x - -v ~ 2 v x , 

3 o 
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and we can reconstruct the corresponding conserved charge using the homotopy formula 



H^= I dx I d\v U -^- AV] 

ov 



to finally obtain 



fl-M = J dx (~\v\l x + ^v 2 vl - 2vv 2 x vl - Iv-'v^j , (48) 

which is also if we make the substitution (|2ip to the right hand side of (|48p . Computationally, 
this is easier than to directly perform the calculations in the u variable. 

From (|47p and (|48p we obtain the following second Hamiltonian structure for the OVF equation 

an 

OU 

V 2 = VD- 2 v- 2 D x v- x D-*v- 1 D x v- 2 D- 2 , 

tt it ( 3 2 ^23 ^22 —1 fi i 

^ = "2 o -«3x- g« v 2x + 4vv x v 2x + -V v x J , 



where u is given by (|2ip . However, the second Hamiltonian formulation for the OV equation (|2|) 
can be written formally as 

u t = V 2 -^- , V 2 = 9D x 2 v~ 2 D x v- 1 D x 3 v- 1 D x v- 2 D x 2 , H ± = ~\f dxu ™ ' ( 49 ) 

where we have used the trivial conserved charge (fT9l) with —IbH^jbu = D 2 • 1 and D~ 2 D 2 -1 = 1. 
This same charge and behavior appears in the pull-back of the Harry Dym equation to the Hunter- 
Saxton equation [27] and is due to the nonlocality of the Hamiltonian operator T> 2 . 



From the CDGSK conserved charge = J dy z we get 

dxv~ l v 2 , 

for the OV equation. In fact, from (|18p and ([49 p we have the recursion operator 

R = V 2 V^ = QD^v-^DbV^D^v^DxV^D- 3 , 

and this could be the starting point to generate a hierarchy of equations and charges for the OV 
system of equations which we will explore in a future publication. 
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